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Abstract. We compare the results of three binary population synthesis codes to un- 
derstand the differences in their results. As a first result we find that when equalizing the 
assumptions the results are similar. The main differences arise from deviating physical 
input. 



Introduction. Binary population synthesis (BPS) codes enable simulations of the en- 
tire evolution of a large number of binary systems. BPS codes are aimed for the study 
of the formation and evolution of astronomical populations, e.g. novae, X-ray binaries 
and type la supernovae (SNela). However, considerable differences arise in the predic- 
tions of e.g. the SNela rate. In this research we compare the results of three different 
BPS codes and investigate the importance of the different assumptions in these codes. 



Method. To compare the results of various BPS codes, we look at two populations 
at their time of formation: Single WD systems with a non-degenerate companion and 
double WD systems. We analyze the evolution paths towards the population of binary 
WDs, in order to have a complete -within the limits of this project- overview of the 
similarities and causes of the differences in the results of the BPS codes. As a first test 
we equalize the assumptions. These assumptions are not necessarily considered to be 
realistic. We simulate a grid with the same initial distributions (in primary mass, mass 
ratio, separation). We assume conservative mass transfer to all types of stars and for 
the common envelope evolution we use the prescription based on Webbink ( 1984]). We 
do not consider magnetic braking, tides and wind accretion. For this paper we only 
consider the formation of WDs with a mass higher than 0.48 M©. Three different BPS 
codes are compared: 



Binary _c: Based on |Hurley et al.| ( [20001 [2002), with updates described in |Izzard 
et al.|(|2006|) and Claeys et al. (in prep.) 



SeBa: Based on Portegie s Zwart & V erbunt (1996), with updates described in 
Nelemans etaLlfl2001D and |Toonen et al.| ( |2012D " 



Brussels code: Based on 



scribed in Mennekens et a 



De Ponder & Vanbeveren| ( |2004[ ), with updates de- 
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Figure 1. Mass of the primary WD versus the separation at the moment of for- 
mation of the second WD, given by the different BPS codes. The numbers of the 
systems are normalized to the total number of the systems forming 2 WDs. 



First results. 

• Similarities: The results of single WDs with a non-degenerate companion are 
very similar. At time of formation of the WD the same regions in separation and 
mass can be distinguished. In the population of double WDs three main regions 
can be distinguished and are similar in the results of the different groups (see 
Fig. [T]). A region with log a > 3, which are the non-interacting systems. A 
region with log a < 1.5, which are systems resulting from a first mass transfer 
phase, which can be stable or unstable, and a second phase which is unstable. In 
the results of Binary _c and SeBa there is a region around log a equal to 1, which 
are systems resulting from a first phase of unstable mass transfer and a second 
phase which is stable. The percentage of double WD systems compared to the 
simulated systems is similar for the different groups, namely 9% for Binary _c, 
for SeBa 10% and 14% for the Brussels code. 

• Differences: The predicted population of binary WDs are similar (see Fig. [I]), 
however small differences in mass and separation can be noticed. These are 
due to variations in the input physics: initial-final mass relation, the criterion to 
determine whether mass transfer continues stable or unstable, the exact duration 
and rate of mass transfer, the rate and angular momentum loss of the wind and 
the prescriptions to describe He-star evolution. 

Conclusion. Various BPS codes give similar results for the formation of binary WDs, 
when the assumptions are equalized. The differences in their results are not due to 
numerical differences, but can be explained by deviations in the input physics. 
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